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Abstract: We study the topological charge fluctuations of an SU(2) lattice gauge theory 
containing both Nf = 2 and 4 flavors of Wilson fermion, at low temperature with non- 
zero chemical potential fi. The topological susceptibility, xt, is used to characterize 
^ ' differing physical regimes as fi is varied between the onset of matter at /io and and color 
■ deconflnement at fid- Suppression of instantons by matter via Debye screening is also 
O^. investigated, revealing effects not captured by perturbative predictions. In particular, 
! the breaking of scale invariance leads to the mean instanton size p becoming /x-dependent 
I in the regime between onset and deconflnement, with a scaling p oc over the range 
fio < < fJ'd, resulting in an enhancement of xt immediately above onset. 

1 Introduction 

Lattice studies of matter at non-zero baryon density are hampered by the 'sign prob- 
lem', which arises when a quark chemical potential term /i is included in the Euclidean 
^ ■ QCD action. The resulting complex nature of the fermion determinant precludes a pos- 
^ ■ itive deflnite probability measure and computational techniques based on importance 
. sampling break down. A gauge theory which is accessible to Monte Carlo simulations is 
^ ' QC2D, based on gauge group SU(2), describing "two color matter". In QC2D, quarks be- 
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long in the pseudoreal 2 representation of SU(2) which can guarantee a positive deflnite 
measure. 

Studies of two color matter have been performed utilising a number of fermion for- 
mulations. The series of works obtained from simulations involving two and four flavors 
^ ! of Wilson fermion [H I21 E] have revealed a scenario in which, as fj, is increased, baryonic 
matter forms at an onset fio = ^■n/'^ whereupon the matter then exists in a superfluid 
state with a progression from a dilute gas of tightly-bound diquark pairs to degenerate 
quark matter, culminating in color deconflnement at around /i ~ l.lm^. This Letter 
supplements this picture with an investigation of topological effects observed on the 
same lattice conflgurations. 

The topological charge density qt may be deflned in terms of the Yang-Mills held 
tensor as ^ 

(It = i^^F^^F^y (1) 



with F^y = \e^ypaFpfj. The action is minimised when the condition = ±F^i/ is satis- 
fled. The observable measured to study topological charge fluctuations is the topological 



1 



susceptibility, xt, defined as 

XT = (2) 

where Q = J d'^xqr and V = J d^x. Using large-A^c methods xt is estimated by means 
of the Witten-Veneziano formula |H |5] 

XT = ^(m^ + m^;,-2my (3) 

to be (180 MeV)^ in the SU(3) gauge vacuum. Simulations of hot two color matter with 
two flavors of staggered quark (equivalent to Nf = 8 continuum quark flavors) have 
shown this quantity drops sharply at the deconfining temperature and have suggested 
this also happens at non-zero chemical potential [SI [!]■ When xt is measured as a 
function of a/i, the susceptibility remains constant before dropping dramatically at a 
critical chemical potential corresponding to both deconfinement and chiral symmetry 
restoration. 

In a semi-classical picture toplological charge is localised on four-dimensional objects 
called instantons, which are solutions of the sef-dual condition for a local minimum of 
the action [8]. Another observable of interest is the size of an instanton p. This is a 
measure of the extent to which the gauge field action is localised. For classical Yang- 
Mills instantons the size may be considered arbitrary due to scale invariance and so p 
does not depend upon the action, and vice versa. However, in the quantum vacuum 
scale invariance is broken, and the typical size of an instanton is estimated to be in the 
region of 0.3 fm P [TO] . 

In dense matter, Debye screening of color charge leads to instanton suppression [12] . 
Perturbative calculations [13] predict that instanton number at large chemical potential 
should go like 

n{p) = n{p = 0) exp (-Nf{ppf) . (4) 

Therefore, as the number of quark flavors Nj is increased, instantons should be sup- 
pressed and Xt should decrease. It should also be expected that, if the average instan- 
ton size p is indeed fixed, then the extra matter present as p is increased will screen the 
topological charge and suppress xt still further. 



2 Methodology 

In order to explore instanton effects on a lattice we replace the continuum topological 
charge density qx ( H]) with its lattice counterpart 

(1l{x) = ^;^tf,yp^i:i{U^^{x)Up^{x)) (5) 

where U^^y^x) is the product of link variables around a plaquette at site x in the p — u 
plane [H]. The charge density is thus measured by taking the trace of the product of 
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two orthogonal plaquettes. The total charge Ql is obtained via Ql = J2x1l{x). Within 
each configuration, the peaks due to the presence of instantons (whose structure may 
extend over a scale p ^ a, where a is the lattice spacing) are mutated by short scale 
(0(a)) fiuctuations. Such UV fiuctuations are highly undesirable as they contribute to 
the total charge but obscure the 'real' instantons, and so the measured susceptibility 
can be an overestimate pS]. The lattice topological susceptibility xl = (QV) differs 
from the continuum value by both a multiplicative factor Z and an additive one M: 

XL = Z^a\T + M. (6) 

Z and M depend on several factors including the quark mass, the inverse coupling (3 
and the choice of fermion operator [7]. In general, on the lattice, Z ^ 1 and the charge 
Ql is not integer-valued. The challenge is to minimise the unwanted, short distance 
contributions while in the process recovering the continuum value in an unambiguous 
fashion. 

for a given configuration of gauge fields is calculated by means of Eqn. ([5]). The 
effects of UV fluctuations are minimised by cooling [16], whereby a new conflguration 
is generated from the old by visiting lattice sites in turn and minimising the action 
locally. Repeating this succesively has the effect of smoothing out fluctuations and 
revealing the underlying topological structure in the gauge flelds. By prudent use of 
cooling, the multiplicative factor Z — 1 as the unwanted fluctuations are eliminated. 
However, excessive cooling eliminates not just the UV fluctuations but will also shrink 
and ultimately eradicate the 'real' instantons. If cooling shrinks an instanton until its 
size p < a then it 'falls through' the lattice and some of the topological information is 
lost. If only larger instantons contribute to the total charge then there is a tendency 
to underestimate Qt- Information can also be lost as too much cooling has a tendency 
to annihilate instanton - anti-instanton pairs. The total charge may remain the same 
but the charge density is reduced. Therefore, it is vital that good control of the cooling 
process is maintained. 

The additive constant M may be dealt with by equating it to the value of the 
topological susceptibility in the Qt = sector, setting M = xo = Xt{Q = 0). As we 
have no prior knowledge to suggest that our ensemble is in the trivial sector we must 
modify Eqn. (|2]). In the non-trivial sector M can be eradicated by redeflning 

. (7) 

Thus, by measuring the charges on a number of cooled fleld conflgurations with Z ~ 
1 and calculating xt by means of ([7]), the physical topological susceptibility can be 
extracted from the lattice one. Henceforth, we discard the references to lattice values 
via our subscripts L and merely label x and Q with the subscript T. 

The cooling method employed here uses a computer program to read the gauge fleld 
information from each conflguration and then calculate the total action by summing over 
the plaquettes. In general, this is not the minimum action. A point is then chosen and 
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a link variable U^{x) is selected. There are 6 plaquettes with this link in common. The 
code sums the link products, in the form of unitary matrices which form the 'staples' 
bordering the link f/^(a;), resulting in a 2 x 2 matrix V. The matrix V is non-unitary 
and must be renormalised as = (Det V)~^^'^V. Keeping V fixed, the action is then 
minimised by modifying U^{x). By systematically working through the old configuration 
and updating all links f/^(a;) a new configuration is produced with a lower action than 
the original one. This completes the first cooling sweep. By predetermining the number 
of sweeps to be performed, the process repeats automatically and the configuration is 
cooled to the required extent. When cooling is complete, the code then searches through 
the final configuration to find where the peaks of the action are located and FF at these 
points is recorded. Setting a minimum cutoff for FF allows the code to disregard the 
smallest fluctuations. Imposing a second cutoff for the maximum extent of the gauge 
fields inside an instanton minimizes any finite volume effects associated with excessively 
large instantons. Once the required topological information is extracted from the cooled 
configuration, the program then moves onto the next configuration in the ensemble and 
repeats as necessary. 

To find the total topological charge on each configuration, a second program obtains 
the net value of all the peaks of FF from the output of the first, providing a sequence 
of estimates for the fluctuating variable Qt- The topological susceptibility is estimated 
from this using Eqn.(I71). 

One aspect of topological structure that is worth investigating is the size distribution 
of the instantons. Instanton size may be calculated from the peak value of the topological 
charge density using 

6 . 

This classical approximation works reasonably well for large lattice instantons, but for 
smaller ones whose size is of the order of the lattice spacing, corrections of 0{a?) are 
needed. The necessary correction factors for A^^c = 3 were calculated by Smith and Teper 
by cooling a classical instanton and then parametrising the resulting relationship between 
Q and p The computational method employed in this study involved reading the 
peak values of the charge from the lattice configurations and then applying iterative 
bisection to find a value for p which satisfied Eqn. ([H]) to within a predetermined error 
factor e. 

3 Numerical Results 

Information about the topological structure was extracted using two different gauge 
field ensembles. The first was generated on a 12'^ x 24 lattice with y/aa = 0.415(18) 
(o" is the string tension) using Nf = 2 flavors of Wilson fermion at an inverse coupling 
P = 1.9 [2j. The fermion action included a diquark source term aj = 0.04 and the 
value of the hopping parameter k, = 0.168. The second ensemble was generated on the 
same system size, /3 and j using Nf = 4, resulting in a signiflcantly flner lattice with 
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^Jaa = 0.138(4) [3]. This time k, was chosen to be 0.158; both ensembles therefore had a 
matched pion mass m.,ra = 0.68(1). Although we choose to plot several figures in cutoff 
units fia, the horizontal axis could therefore equally be regarded as being calibrated 
in units of fi/m.„, as noted in [3j. In both cases chemical potential was introduced 
via the standard Hasenfratz-Karsch prescription [18]. The minimum cutoff for FF was 
Qcut = 0.02a~^ whereas any instantons larger than one third of the spatial extent of the 
lattice were ignored. 




a|j. 



Figure 1: Topological susceptibility a^XT versus chemical potential for N f = 2 andNj = 4 
performing 10 cooling steps. 

The topological susceptibility of the Nf = 2 and 4 configurations was measured across 
a range ~ 0.25 < afi < 1.1. Fig. [T]shows the behaviour of xt for both theories. In order 
to verify the validity of our approach to cooling, the same ensemble was submitted to 
both 10 and 20 cooling steps. As the extra cooling had little effect on the signal, we are 
confident we are characterising the underlying topology satisfactorily. 

In the Nf = 2 case the signal remains fairly consistent across the range studied. There 
are possible signs of some minor downward fiuctuations at afi ~ 0.4 and a/i ~ 0.5, but 
beyond a/x = 0.5 the data remains fiat. The Nf = 4 data obtained on a finer lattice are 
more interesting. Fig. [T] suggests that the extra fiavors have suppressed the instantons 
with the peak value of Xt^~'^ ~ ^■^Xt^~'^- Moreover, the suppression of Xt^~ the 
lowest densities is much greater can be explained purely by a change in Nf, and its 
behaviour with increasing fi is not as expected if the relation in Eqn.(jl]) is correct (and 
assuming p to be independent of Nf); x(yu) initially increases instead of being suppressed. 

At larger densities the situation changes. A comparison of Xt^~^ with the Polyakov 
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Figure 2: T/ie suppression of xt coinciding with the rise in (L) for Nf = 4. Note (L) 
has been rescaled for clarity. 



loop (L) from [3] over the range 0.7 < afi < 1.2 is shown in Fig. |2l It illustrates nicely 
how the fall in xt coincides with the rise in (L). The Polyakov loop begins to rise from 
zero at afi ~ 0.8, whereas xt starts to fall noticeably just a little later at afi ~ 0.9. 
Ref. [3] identified a "deconfining" value of chemical potential a fid ~ 0.75 based on the 
behaviour of (L); Fig. [2] suggests deconfinement in dense matter is accompanied by 
suppression of topological fluctuations. 

It is worth noting that although effort has been made to make a direct comparison of 
the Nf = 2 and Nf = 4 data sets, in reality they are distinct theories. When comparing 
both ensembles it is important to remember that the physical volumes (ie. as measured 
in string tension units) of the two lattices differ by a factor of ~ 3^. Hence it is difficult 
to make any useful quantitative comparison of xt for the two theories directly. It is 
more useful to rescale xt in each case as some dimensionless parameter and then make 
a comparison of the two. Fig. [3] shows the topological susceptibility rescaled and plotted 
as the fourth root of a'^XT divided by the square root of the string tension a'^a. From this 
it is possible to compare the results to ([3]), which implies Xr^ = 180 MeV. For Nf = 2, 
xl^^a^/'^ = 0.3493 ± 0.0076, and for Nf = 4, Xr V^^^^^ = 0.4837 ± 0.0198. Assuming 
a = (440MeV)2 leads to 
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156 ± 3MeV Nf = 2; 
^ 213 ± 9MeV Nf = 4. 

Both results are in the range suggested by the Witten-Veneziano formula. That the 
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Figure 3: The dimensionless quantity Xt/'^^ versus ^/cr^ for 2 and 4 flavors. 

value for Nf = 2 is smaller than the Nf = 4 value by about 25% hints that much of the 
topological information is being missed due to the coarseness of the lattice. 

The instanton size distribution was plotted for both the two and four flavor config- 
urations. To compensate for inequalities in the sizes of the data sets, each distribution 
curve is normalized, keeping the area under the curve consistent. Using the same scale 
determination that led to Eqn. (|9]), for Nf = 2 (see Fig. H]) at a/i = 0.25, the majority 
lay within the range 0.18 < p < 0.5fm, with the average size being p ~ 0.28fm. This 
compares well with the phenomeno logically derived value p ~ 0.3 fm [17]. The sharp 
cutoff at p = a is where instantons smaller than this 'fall through' the lattice and do 
not contribute. Similar distribution curves are also plotted for higher values of p to see 
if there is any effect with increasing density. While there is a hint that larger instantons 
are shghtly suppressed at larger p, no significant /i-dependence is observed and all the 
curves are qualitatively the same. 

By contrast. Fig. |5] shows how the instanton size distribution evolves with p for 
Nf = 4. At ap = 0.25 the distribution is fairly uniform. As p increases the number 
of large instantons falls as that of smaller-sized instantons rises, and the distribution 
becomes taller and narrower, with a peak at ap = 0.9 of p ~ 1.4fm. For ap = 1.0 
the curve has rapidly flattened and has a very similar profile to that for ap = 0.45. 
The prevalence of small-sized instantons drives down the average instanton size. The 
cutoff at large p for ap = 0.25, 0.45 results from the constraint on the maximum possible 
instanton size. Such a filter on p is needed to minimize finite volume effects and to 
prevent instantons overlapping one another. As p increases, p decreases and this becomes 
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Figure 4: Normalized instanton size distribution at three different chemical potentials 
for Nf = 2. The average size p 0.3 fm. 

less of an issue but at low /i, where there is a greater number of large-sized instantons, 
it seems likely that some topological information is lost due to the IR cutoff. 

When p is plotted as a function of fi (see Fig. [6]) it is apparent that the iVj = 4 
instantons are shrinking between onset and deconfinement, after which there is a sharp 
rise in size. This rise coincides with the rapid flattening of the a/x = 1.0 distribution curve 
recorded in Fig. [5l In the range /io < < p^d, empirically p oc p,~^ and this behaviour 
is plotted along with the data. This contrasts markedly with the Nf = 2 data, where p 
appears to be almost /^-independent. At most there is a very gentle monotonic fall in p 
for p, > p,o. A consideration of the behaviour of the size distribution in Fig. [5] suggests 
why that might be. The most significant /i-dependent effects are seen for small-sized 
instantons with p < 3a^-f^^ ~ a^f^"^. Thus, the Nf = 2 lattice may well be too coarse 
for this detail to be seen. 

It is also of interest to compare the behaviour of p(/i) for A'^^ = 4 with the predictions 
of thermal field theory. The perturbative result (j4]) of [13] implies that, for constant p n 
is suppressed by a factor e~'^fp ^ as density is increased. Thus, xt should be suppressed 
by increasing /i. In Fig. [T] Xt^~^ rises as p, increases, which seems to be incompatible 
with the perturbative result. However, if we take into account the non-perturbative 
information on p(/i) oc in Eqn.(jl]) we find 
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Figure 5: Behaviour of the normalized Nj = A instanton size distribution at various fi 
values. 

The parameter A includes n{fi = T = 0) plus thermal contributions due to the fact 
that we are working at low but non-zero T. The parameter B includes the constant of 
proportionality for instanton size as a function of fi multiplied by a factor of Nf. When 

a function of this form is plotted along with Xt Fig- d the range iJ,o < fJ. < fJ-d 

we see a fair correspondence between the two. The best fit is found with A = 1.8255 
and B = 0.2231. This suggests that the enhanced topological fluctuations observed 
in baryonic matter at moderate density are a direct result of the /x-dependence of the 
instanton scale size. 

4 Conclusion 

In this Letter we have presented the first exploratory study of topological fluctuations 
of non-abelian gauge fields in cold dense baryonic matter, using ensembles generated for 
a range of with both Nf = 2 and Nf = 4. 

For Nf = 2 the topological susceptibility remained fiat across the whole range stud- 
ied. It is likely that the lattice is too coarse to be able to capture the topological 
detail adequately. The fact that for the two flavor ensemble xt was measured to be 
Xt — (150 MeV)^ suggests that a lot of topological information is falling through the 
lattice and being lost. 

With an increase in the number of flavors to Nf = 4, the resulting finer lattice was 
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Figure 6: Average instanton size versus chemical potential for Nj = 2 and for Nj = 4 
fitted with function f{fi) oc for fio < ft < fid- 



able to expose more detail about the distribution of instanton size and its evolution 
with /i. Once chemical potential is increased beyond onset, the instanton size becomes 
density-dependent. The large instantons found at low /i are suppressed as /i is increased, 
driving down the average size to a minimum of p ~ 0.14 fm at around afi ~ 0.9. The 
smaller instantons result in reduction in screening of topological charge fluctuations, so 
that the topological susceptibility xt initially rises for > Ho- In the deconfined phase 
> fJ-d, however, the average instanton size rises sharply, and xt is suppressed. 
While it is tempting to ascribe the differences observed between Nf = 2 and Nf = 4 
entirely to the different lattice spacings, as measured in string tension units, we should 
remain mindful that they are two different theories; in particular the thermodynamics 
studies of [21 13] reveal that for Nf = 2 the regime just above onset is weakly-interacting 
and dilute, apparently well-described as a non-relativistic Bose gas of tightly-bound 
scalar diquarks. By contrast, matter with Nf = 4 appears relativistic and strongly- 
interacting for all n > Ho- A systematic study of the /i-dependence of topological 
fluctuations in QC2D, therefore, must await the generation of gauge ensembles on a 
finer lattice. 
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